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Abstract 

An infinite square well with a discontinuous step is one of the simplest systems to exhibit non- 
Newtonian ray-splitting periodic orbits in the semiclassical limit. This system is analyzed using 
both time-independent perturbation theory (PT) and periodic-orbit theory and the approximate 
formulas for the energy eigenvalues derived from these two approaches are compared. The periodic 
orbits of the system can be divided into classes according to how many times they reflect from the 
potential step. Different classes of orbits contribute to different orders of PT. The dominant term 
in the second-order PT correction is due to non-Newtonian orbits that reflect from the step exactly 
once. In the limit in which PT converges the periodic-orbit theory results agree with those of PT, 
but outside of this limit the periodic-orbit theory gives much more accurate results for energies 
above the potential step. 
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I. INTRODUTION 



Periodic-orbit theory is one of the most interesting developments in the study of the 
relationship between classical and quantum mechanics [T]. The centerpiece of periodic- 
orbit theory is the Gutzwiller trace formula which relates the quantum density of states 
to properties of the periodic orbits in the classical system. In general the trace formula 
gives only approximate results, becoming exact in the limit h — > 0. However, for some 
systems the periodic-orbit theory is exact even for finite H [2]. Periodic-orbit theory has 
found particularly fruitful application in the study of quantum systems with chaotic classical 
counterpart [TJ EH5] . 

Recently periodic-orbit theory has been extended to the case of ray-splitting systems [6] . 
In most systems the wavelength of a quantum particle in the semiclassical limit (h — > 0) is 
small compared to all relevant length scales in the classical system. In this case the quantum 
wave equations reduce to ray equations and the particle obeys Newtonian mechanics. In 
ray-splitting systems, however, the potential changes significantly even on length scales 
that are small compared to the wavelength of the quantum particle in the semiclassical 
limit. This will occur if there is a discontinuous change in the potential within the region 
accessible to the particle. In the semiclassical limit of a ray-splitting system a particle 
will follow Newtonian mechanics everywhere except at the discontinuous boundary. At the 
discontinuity the particle may be transmitted across the boundary or reflected from it, so 
the ray splits into two parts. The non-Newtonian periodic orbits that result from this 
ray-splitting can influence the quantum dynamics. Recent computational and experimental 
studies of ray-splitting systems have revealed signatures of non-Newtonian periodic orbits 
in the Fourier transform of the quantum density of states [71412]. the distribution of level 
spacings [8j [131 E]> an d the scarring of energy eigenstates [TJ El [101 IT2] . 

Periodic-orbit theory is typically applied to systems with chaotic classical dynamics. The 
quantum versions of these systems are generally not amenable to exact solution, or even 
approximation methods like perturbation theory, and must be studied numerically. Ray- 
splitting systems are interesting in this regard because their classical dynamics may exhibit 
some properties of chaos (like an exponential proliferation of periodic orbits with increasing 
period) but perturbation theory may still be applied to the quantum dynamics of the system 
(at least for certain parameter regimes). For some of these systems exact formulas have been 
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found which give the quantum energy eigenstates as a sum over the classical (Newtonian 
and non-Newtonian) periodic orbits. These cases allow for a direct analytical comparison 
between the results of perturbation theory and those of periodic-orbit theory. The goal of 
this study is to carry out this comparison for the simplest possible system: the asymmetric 
infinite square well (AISW). 

The AISW consists of an infinite square well of width 2a with a discontinuous step of 
height Vq at the center of the well [T5HPT] . The potential energy function is 



f 

oo, Ixl > a 



V{x) 



0, -a < x < (1) 
Vq, < x < a. 



A plane wave with energy E > Vq incident on the boundary at x = may be reflected with 
probability r 2 or transmitted with probability 1 — r 2 , where 



i - yi - V./E 



i + v 1 - V »I E 

Note that r does not depend on % so these non-Newtonian reflections persist in the semi- 
classical limit and the classical dynamics contains non-Newtonian periodic orbits such as 
those shown in Figure [TJ The orbit iV is the Newtonian orbit that moves back and forth 
between the hard walls at x — ±a. The orbit L is a non-Newtonian orbit that reflects when 
it reaches the boundary at x = from the left, so this orbit is confined to the left side of 
the well. The orbit R is confined to the right side of the well, reflecting when it reaches the 
boundary at x — 0. These three basic orbits can be combined to form an infinite variety of 
other periodic orbits. Non- Newtonian orbits are also possible for E < Vo, but in this case 
the orbits are "ghost orbits" that explore the classically-forbidden right side of the well. For 
perturbation theory to be valid the energy of the particle must be greater than Vq, so there 
will be no need to consider the case E < Vq in this paper. 

Section presents the results of a second-order (Rayleigh-Schrodinger) time-independent 



perturbation theory (PT) analysis of the AISW. Section III details the application of 
periodic-orbit theory to this system, with the goal of providing an approximate formula for 
the energies that is comparable in accuracy to the second-order PT results. The applica- 
tion of periodic-orbit theory to the AISW clearly shows that periodic orbits with different 
numbers of reflections at x = contribute to different orders of PT. The second-order PT 
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correction is due predominantly to non-Newtonian orbits that have a single reflection from 



the potential step. Section IV presents a comparison of the results of PT and periodic-orbit 
theory to the exact energy eigenvalues. The approximation derived from periodic-orbit 
theory is shown to be more accurate than that derived from PT. Section [V] provides a 
summary and discussion of the results. Most of the details for the PT and periodic-orbit 
theory calculations are given in Appendices at the end of the paper. 

II. PERTURBATION THEORY 

Standard (Rayleigh-Schrodinger) time- independent perturbation theory proceeds by writ- 
ing the full Hamiltonian for the system as H — if + H' where Ho is the Hamiltonian of the 
"unperturbed" system and H' is the perturbation. For the AISW the unperturbed system 
is simply an infinite square well with potential function 

. , oo, \x\ > a 
V(x) = { (3) 

I 0, |x| < a 

where the width of the well is 2a. This system is discussed in almost any textbook on 
elementary quantum mechanics [IE]. The energy eigenvalues of this system are 
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E ( ) = vTin 
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where m is the mass of the particle in the well. The wave function for the energy eigenstates 
are 

(o)/„\ _ 1 ( mr(x + a) 



The perturbation for the AISW is then 



H' 



0, x < 

Vb, < x < a (6) 
0, x > a 



where Vq is the height of the potential step inside the well. The infinite square well has no 
degeneracies, so non- degenerate PT can be applied to find approximations for the energy 
eigenvalues of the AISW. The first order correction is 

= MW^) = Voj o a s^ (^|^) dx = \ (7) 



This correction simply shifts each energy eigenvalue upward by half the height of the poten- 
tial step. The second order correction is given by 
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The second order correction is evaluated in App. [Aland the result is 
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where 

3, n is even 
—1, n is odd 

and 0(g(n)) indicates a function f(n) that is less than some constant times g(n) for all n 
greater than some value n = C. So to second order the energies of the AISW are given by 
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The 0(1 /n 3 ) terms can be ignored for sufficiently large n. 

To simplify the notation it is helpful to introduce a dimensionless constant 
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Note that this constant provides information about the size of the potential step Vq relative 
to the ground state energy of the infinite square well E[°\ since 
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In terms of this dimensionless constant the AISW energies are approximately given by 
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Before any attempt is made to use the results of PT one must carefully consider whether 
or not the perturbation expansion will converge. The requirement for the rapid convergence 
of the perturbation expansion is [19] 
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It is clear that the left side of Eq. 15 will be largest when k = n ± 1, in which case 

8a 2n ± 1 =f (-l) n 
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So the condition for the convergence of PT reduces to 



« 1. (17) 
n 



III. PERIODIC ORBIT THEORY 

Dabaghian and Jensen derived an exact formula for the energy eigenvalues for the AISW 
in terms of an infinite sum over the periodic orbits of the classical system |T7j. The procedure 
they used does not converge for E < Vq [20], but a convergent semiclassical formula for 
E < Vo can be obtained by accounting for ghost orbits (orbits that exist in the right side of 
the well, which is forbidden in Newtonian mechanics) [2, 20J. The analysis in this paper is 
limited to large n for which E n > Vq, so the effect of ghost orbits will be ignored and the 
simpler formula of Ref. [T7| will be used. 

Dabghian and Jensen give their semiclassical formula in terms of the discrete quantum 
actions S n . These quantum actions are related to the energy eigenvales of the system via 
the definition of the classical action length for a particle with energy E moving across the 
well: 

S(E) = aV2mE + a^2m(E - V Q ). (18) 
This equation can be inverted to find the energy as a function of action 



E(S) = < ffl 8 + y (19) 



and the quantum energy eigenvalues are then given by E(S n ). 

The semiclassical formula for the reduced quantum action (s n = S n /K) associated with 
the nth energy eigenstate is [T7] : 

-7r(n+l/2) i /■ 7 r(n+l/2) fiu 
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where the variable s is the reduced classical action length s = S/%. N(s) is the Weyl average 
for the density of states with E > V in the AISW, given by 

= S -~\. (21) 

7T I 



The first integral in Eq. 20 can be evaluated and the formula for the reduced actions can 
be simplified to 

s n = n(n - LO n ) (22) 

where 

I _ Mn+l/2) J±v 

w n E-ImV -JL e *»»Pd8. (23) 
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The index p labels all of the fundamental periodic orbits of the system (orbits which are not 
simply repetitions of other periodic orbits). The index v accounts for repetitions of these 
orbits. The classical action of each periodic orbits is given by 

S p = n L 2a\ / 2mE + n R 2a^2m(E - V ) (24) 

where til and are the number of times the orbits passes back and forth across the left 



and right sides of the well, respectively. The quantity s p which appears in Eq. 23 is then 
the reduced classical action for a periodic orbit: 

S p ( . 2a(n L - n R ) 
s p = — = {n L + n R )s^ . (25) 

The factor A p is a weighting factor given by 

A p = (_i)x(p) r «^(p) t r(p) (26) 

where a(p) counts the number of times the orbit p reflects from the barrier at x — 0, r(jp) 
counts the number of times the orbit transmits through the barrier, and x(p) counts the 
combined number of reflections from the hard walls and right reflections from the boundary 
(each of which results in a sign change in the particle's wave function). The reflection 
coefficient r is found by substituting Eq. [19] into Eq. [2] to find 
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with a defined as in Eq. [T2j The transmission cofficient t is given by 
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t 2 = \-r 2 = \ -. (28) 



A more detailed discussion of this formula can be found in Ref. [T71 



Substituting Eq. 22 into Eq. 19 provides a formula for the energy eigenvalues of the 
AISW in terms of n and u;„: 
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A comparison of this result with the PT formula (Eq. 14) shows that the zeroth and first 



order PT terms are also present in Eq. 29 These terms arise from the non-oscillatory 



part of Eq. 20 which is associated with the Weyl average (Eq. |2lj). To compare the 



periodic-orbit theory result to the second-order PT correction it is necessary to evaluate the 



oscillatory term u n using Eq. 23 This formula involves a sum over all (Newtonian and 



non-Newtonian) periodic orbits of the classical system, including all repetitions of periodic 
orbits. The Newtonian orbit and its repetitions will be addressed first. 

A. Newtonian orbits 



The Newtonian periodic orbit passes back and forth across the entire well from x = —a 
to x = a without reflecting from the potential discontinuity at a; = 0. For this orbit 
n>L = n-R = 1, and Eq. 25 then gives s p = 2s. Since this orbit reflects once off of each of 



the two hard walls, but does not reflect from the potential step, it is clear that a(p) = 0, 
r(p) = 2, and xip) = 2- If the contribution of the Newtonian orbits and its repetitions to 
uj n is designated by u nt0 then 
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The details of this calculation are given in App. B 1 Now u n appears twice in Eq. p29l once 



in a term involving nuj n and again in a term involving u 2 . From Eq. 30 it is clear that 
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(31) 
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Therefore the Newtonian orbit and its repetitions make a contribution to the energy formula 
in Eq. 29 that is smaller than other terms that have been ignored in that equation. At the 
level of approximation given by Eq. 29 the Newtonian orbit can be ignored. 



B. Non-Newtonian orbits 



Since the Newtonian periodic orbit makes no significant contribution to the energy for- 



mula in Eq. [29] it is necessary to examine the contributions of non-Newtonian periodic 
orbits. It is convenient to divide the non-Newtonian orbits into classes based on how many 
times each orbits reflects from the potential step at x — 0. If pk designates a periodic orbit 
with k reflections at the step then cr(p k ) = k, so the contribution to u n from orbits with k 
reflections is 

»7r(n+l/2) f2 a \ k ( 4 a 2\ 
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where the sum is taken over all orbits with k reflections at the step. Note the absence 



of the sum over u, which was used to account for repetitions of orbits in Eq. [23] When 
calculating the contribution of the orbits pk to u n it is unnecessary to consider repetitions 
of these orbits, since the repetition of an orbit with k reflections would be an orbit with 2k 
reflections which would belong to a different class. 

Examination of Eq. [33] reveals the way in which different classes of orbits contribute to 
different orders of perturbation theory. The formula for contains a term proportional 
to a k and other terms that involve higher powers of a. The feth order correction from 
perturbation theory is always proportional to a k so periodic orbits with k reflections can 
only contribute to perturbation theory corrections of fcth order and higher. The Newtonian 
orbit (and its repetitions) can contribute to all orders, but we have seen above that the 
zeroth and first order corrections (as well as part of the second order correction) come from 
the Weyl average while the contribution of the Newtonian orbit to the oscillatory term uo n 
has a negligible effect on the second order term. Single reflection orbits can contribute, in 
principal, to first order and higher corrections. Orbits with three or more reflections cannot 
contribute to the second order correction from perturbation theory, so they need not be 



considered in comparing Eq. 29 to Eq. 14 



A closer look at Eq. [33] reveals that two-reflection orbits make only a negligible con- 
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tribution to the second order PT correction. The variable of integration, s, in Eq. 33 
approximately equal to ttti throughout the range of integration, so 
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So c<j n 2 G 0(l/n 4 ) and the contribution to the energy formula from two-reflection orbits is 
0(l/n 3 ). Since Eq. 14 already ignores terms of this order the contribution of two-reflection 
orbits can also be ignored. 

These arguments imply that the dominant term from the second order PT correction 
must come from the contribution of single reflection non-Newtonian periodic orbits (in com- 



bination with the 4a 2 / (n A n 2 ) term in Eq. 29 which comes from the non-oscillatory part of 



Eq. 22). This result is shown explicitly in the next section. 



C. Single reflection orbits 



Appendix B2 provides a detailed calculation of u! n ,i- The result is 
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Since the contribution of the Newtonian orbit and its repetitions is u) n ^ G 0(l/n 5 ) and the 
contribution of ^-reflection orbits is oj n ^ G 0(1 /n 2k ) as shown above, then 
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So in comparing Eq. 29 with Eq. [L4] only the contributions to uo n from single reflection 



non-Newtonian orbits need to be considered. Inserting Eq. p5Hn place of co n in Eq. 29 gives 
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The perturbation theory approximation of Eq. 14 is only valid if oc/n « 1. In this case 



the sine function in Eq. 37 can be approximated as 
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Inserting this result into Eq. 37 and combining like terms gives 
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with 7 defined as in Eq. 10 So the formula for the approximate energies derived from 
periodic-orbit theory, including only the contributions to u n from single reflection orbits, 



matches the formula derived from second order PT (Eq. 14) in the regime in which PT is 
valid. 



IV. COMPARISON OF APPROXIMATIONS AND EXACT ENERGIES 



Although the approximations given in Eqs. 14 and 37 are equivalent in the limit a/n « 



1, it is instructive to examine how well each of these approximations performs outside of 
this limit by comparing these approximations to the exact energy eigenvalues of the AISW. 
For E > Vq the energy eigenstates for the AISW are of the form [15] 



, Asm\Q(x + a)} for — a < x < 

V(aOH ~ (40) 

B sin [q{x — a)] for < x < a, 

where Q = \ / 2mE/h and q = y^2m(E — Vq)/Tl. Requiring ip and dijj/dx to be continuous 
at x = leads to the energy eigenvalue equation: 

Q cos(Qa) sin(g6) + q cos(gfr) sin(Qa) = 0. (41) 

The nth energy eigenvalue will lie in the open interval (E n , E n+ i) where 

E n = E(S)\s=hw(n-l/2) (42) 



with E(S) given in Eq. 19 [17J. A simple bisection method can be used to rapidly solve Eq. 
41 on this open interval to find the value for E n [21J. This procedure can be automated to 
find any number of eigenvalues. 

Once the eigenvales are calculated they can be compared with the approximations from 
periodic-orbit theory and PT. Figure [2] shows the numerically computed energy eigenvalues 



for the AISW and the approximate values derived from second order PT (Eq. 14 ) and from 



periodic-orbit theory using only single reflection orbits (Eq. 37). The parameter values for 
the data shown are (in scaled units): a = 3, Vq = 100, m = 1/2, and H = 1 so that the 
dimensionless quantity a = 450. Thus for all of the data shown in Fig. [2] a/n > 20 and 
we might not expect either approximation to work well. For the parameter values used in 
Fig. [2] there are 10 eigenstates with E n < Vq. Note that neither PT nor periodic-orbit 
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theory provide good approximations for n < 10. The periodic-orbit theory approximation, 
though, is quite accurate for n > 10. The agreement for n = 10 is somewhat surprising 
since the n = 10 state has an energy below the step in the regime where ghost orbits must 
be considered and the version of periodic-orbit theory presented above is not strictly valid. 
For eigenvalues greater than Vq the periodic-orbit approximation is clearly superior to the 
PT approximation, even though the two approximations agree in the limit a/n << 1. This 
general pattern holds true for other sets of parameter values with a > 1. 



V. CONCLUSION 



The asymmetric infinite square well, which consists of an infinite square well with a dis- 
continuous potential step at the center, is an unusual system because its energy eigenvalues 
can be approximated using both perturbation theory (which is usually applied only to sys- 
tems with near-integrable classical dynamics) and periodic-orbit theory (which is usually 
applied only to classically chaotic systems). A comparison of these two approximations re- 
veals that different classes of periodic orbits contribute to different orders of perturbation 
theory (PT). The zeroth and first order terms in the PT expansion, as well as part of the 
second order term, can be derived using only the non-oscillatory term (the Weyl average) 
from periodic-orbit theory. A close examination of the oscillatory part of the periodic-orbit 
theory formula reveals that orbits with k reflections from the potential step can contribute 
only to fcth order and higher terms in the PT expansion. Aside from the terms already 
accounted for by the Weyl average, the Newtonian orbit and its repetitions (which have no 
reflections) make no contribution to the zeroth or first order terms and their contribution to 
the second order term is negligible. Periodic orbits with a single reflection from the potential 
step contribute the dominant term in the second order PT correction. 

Much of the behavior described above is likely to carry over to a wide variety of other 
ray-splitting systems. At high energies the coefficient r for reflection from the potential step 
will be small. Because every reflection contributes a factor of r to the weighting factor A p in 



the periodic orbit sum (Eq. 20), orbits with fewer reflections will generally make larger con- 
tributions to the sum. Orbits with more reflections will contribute only to higher order terms 
in the PT expansion. It also seems that for any ray-splitting system the Newtonian orbit 
will contribute only to second-order and higher terms in the PT expansion. This indicates 
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that single-reflection orbits play a particularly important role in providing a semiclassical 
explanation of the eigenvalue spectrum in ray-splitting systems. 

It is also interesting to note that for the AISW periodic-orbit theory (using oscillatory 
contributions only from single reflection orbits) provides a much more accurate approxi- 
mation to the energy eigenvalues with E > Vq than does standard perturbation theory. 
This raises the possibility of using periodic-orbit theory to find accurate approximations for 
energy eigenvalues in other ray-splitting systems. It would also be interesting to compare 
the periodic-orbit theory approximation with results from other forms of time-independent 
perturbation theory like the Dalgarno-Lewis method [22] or logarithmic perturbation theory 
[23] . High-accuracy analytical approximations for the energy eigenvalues of ray-splitting 
systems could be useful for a variety of applications, such as studying wave packet revivals 
in these systems. 



Appendix A: Second Order Perturbation Theory 

The second order correction from PT is given by Eq. [8] where the matrix element in the 
denominator is 
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The absolute value of this matrix element simplifies to 
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So the second order PT correction (Eq. |8| is 
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The infinite sums in Eq. |A3| can be approximated by first noting that only values of k close 



to n (but differing from n by an odd number) will make significant contributions to these 
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sums. If k = n + 2i — 1 then the first sum in Eq. A3 can be written 
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If n is large then the 0(l/n 5 ) terms can be ignored. Now 
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absolute value. The other sum in Eq. A4 



since each positive term in the sum is cancelled by a corresponding negative term of equal 

is 
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Inserting the results from Eq. |A5| and A6 into Eq. A4 gives 
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Inserting this result into Eq. A3 shows that for even n 
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The sum that appears in Eq. A3 for odd n can be approximated (using k = n + 2i — 1) 
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The first sum on the right hand side of Eq. A9 is identical to the sum in Eq. A4 and thus 



it evaluates to the result given in Eq. A7 The second sum on the right hand side of Eq. 



A9 can be expanded in powers of 1/n and if only the lowest order term is kept the result is 
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Inserting this result into Eq. A3 shows that for odd 
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Combining the results for even and odd n provides the result given in Eq. M above. 



Appendix B: Periodic- Orbit Theory Approximations 



Newtonian Orbits 



The contribution of the Newtonian orbit and its repetitions to u n is given by 

r-7r(n+l/2) / 4^2 \ v 



n Mn-l 

1 x°° -v 1 f f 7T ( n+ 



z i2us ds 



7r(n+l/2) 



e i2us ds - Ava' 



7r(n+l/2) Alvs 



-ds 



O 



7r(n-l/2) 



rr 



The first integral in the right hand side of Eq. Bl is easily evaluated: 

»7r(n+l/2) 



3 l2uS ds 



7r(n-l/2) 



0. 



(Bl) 



(B2) 



The second integral in the right hand side of Eq. Bl can be approximated by changing the 



variable of integration x = s — mr and expanding the integrand in powers of 1/n to find 

t/2 gj2i/(ri7r+a;) 



7r(n+l/2) e i2^s 



7r(n-l/2) 



tt/2 ("T + ^) 4 
1 /-tt/2 



dx 



>*<>■ 1 1 l r V/,. . o ( 



tt/2 



Till J 



n 



(B3) 



The integral on the right hand side of Eq. B3 can then be split into two parts. The first part 



evaluates to zero, since it is equivalent to the integral in Eq. B2 with n = 0. The second 
part involves the integral 



n/2 



xe uvx dx 



tt/2 



2u 



,, v sin(i/7r) — txv cos(wn)) 



ztt(-I) 



v+l 



2v 



Substituting these results back into Eq. |B1| gives 



7T 6 n 5 



y-(-l) 1 
^ v 



O 



TV 



The infinite sum in Eq. |B5| can be evaluated: 

- (-1)^+1 

f=i 



log(2) 



(B4) 



(B5) 



(B6) 



and substituting this result into Eq. B5 gives the result shown in Eq. M0\ above. 
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2. Single Reflection Orbits 



Periodic orbits in the AISW can be represented as sequences of the letters L and R, 
where L represents a back-and-forth motion across the left side of the well and R indicates 
back-and-forth motion across the right side of the well (Fig. [I]). The Newtonian orbit is 
represented by LR (or RL) since it involves motion on one side of the well immediately 
followed by motion on the other side. Examples of orbits with a single reflection are: L, R, 
LRL, RRL, LRLLRLR, etc. The cyclic permutation of the symbols in the representation 
of a periodic orbit just produces a different representation of the same orbit (but starting at 
a different point in the cycle). Therefore all cyclic permutations of symbols are considered 
equivalent. For example, LRL and RLL are really the same orbit. Each group of equivalent 
symbol sequences defines a necklace (a name chosen to conjure the image of arranging the 
symbols in a circle). A particular sequence can be chosen to represent each necklace. 

The symbolic necklaces for periodic orbits with exactly one reflection fall into two groups. 
The first group consists of necklaces that can be represented by some number of repetitions 
of the sequence LR followed by a single L. Table |T] lists several of the necklaces in this group 
along with the values of til + n Ri T {p), an d x(p) f° r each necklace. The bottom row of the 
table provides general formulas for the jth necklace in the group. The contribution of this 
group to u n7 i is designated u n ^L and is given by 

?n «2L Mn+l/2) i / a 2 

~ /! ' V(-l) 2 ^ 1 / e i[(^l)s+2 a /s] ds 
j=l A(n-l/2) s V s J 

9rv Mn+1/2) i(2j-l)s p i2a/s / i \ 

S5> /, / d3+0 {^)- (B7 > 

71 j=l Jw(n-l/2) s \ n / 



Un,lL — 9 



This integral can be written in terms of the variable x = s — mr to give 

r> 00 rn/2 -i 

-Ot . ,,,..-/ I 

Note that 

e i(2j-l)nn = ^_^n ^ 



"n*L - ~ Elm C *W-D- —L^-xy.^te + o(\). (B8) 
n 2 ^ J_ n/2 [n-K + xy \n b ) 



and two of the factors in the integrand of Eq. B8 can be expanded in powers of 1/n to give 



e £f, = e^/(-) ( 1 - ^ ) + O ( ^ ) (BIO) 



n 2 -K 2 J \n 3 
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and 



(nir + x) 2 n 2 7T 2 



1-^+0 
mx ) 



n' 



(Bll) 



Substituting these results into Eq. B8 gives 



n,iL = (-l)" +1 4^Im e* 2 ^V 

j'=i 



00 /-7T/2 

-tt/2 



(2nn + ?2a)x 



n 2 n 2 



e l{2] - 1)x dx + O 



(B12) 



The integral in Eq. B12 can be evaluated to give 



tt/2 

-71-/2 



n 2 ir 2 



c dx 



tt/2 



7T/2 



J(2j-l)x dx 



2wr + i2a f n/2 



n 2 n 2 



xe 



tt/2 



2 cos(j7r) z2n7r — 2a 2 cos(j7r) + (2j — l)7rsin(j7r) 



n 2 7T 2 



+ 



1/ 



2j - 1 wr(2j - 1)< 



(2j - I) 2 

o- 1 



(B13) 



since cos(j7r) = (— l)-' and sin(j7r) = 0. The exponential factor in front of the sum in Eq. 



B12 can be expanded using Euler's formula: 

i2a/(niz) 



/2a \ . . 
cos — + 1 sin 



ymr J 



2a\ 
nn J 



(B14) 



Substituting Eqs. B13 and B14 into Eq. B12 gives 



2a(-l) n+1 



n 2 n 4 



2 sin 



(2a\ 

cos — 



mr 



E 



'_iy+i 



Note that 



and 



E 



' 2? - 1 

l)i+i 



7T 



E 



G 



(B15) 
(B16) 

(B17) 



where G ~ 0.916 is Catalan's constant. Inserting the results for these sums into Eq. B15 
gives 



2a(-l) 



n+l 



n% 4 



7T /2a 
-sin — 
2 V rwr 



4G 

H cos 

nil 



2a\ 
nn J 



O 



(B18) 



The necklaces in the second group of single reflection orbits can be represented by rep- 
etitions of the sequence RL followed by a single R. Table |TT] provides values of n L + ur, 
t(p), and x{p) f° r several of these necklaces as well as general formulas for the jth necklace 
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in this group. The contribution to u> n> i from this group of orbits is designated 0J n> iR and is 
given by 



UJ 



7T Z *■ — ' 

J'=l 



7r(n+l/2) j 



7r(n-l/2) 



1 - 



4a' 



2\ i-i 



3 i[(2j-l)s-2a/a] ds _ 



(B19) 



Note that Eq. B19 is identical to Eq. B7 except for the change a — > —a. So cj n i# can be 



evaluated by simply changing a — > — a in Eq. B18 The result is 



2a(-l) 



n+l 



n 2 7r 4 



7T 



Sill 



2a \ 4G 

— cos 

mi J nn 



2a \ 
nn J 



O 



n' 



(B20) 



The contribution to u n from all of the single reflection orbits is simply the sum of the 



contributions from the two groups: u n i = ^n,iL + ^n,iR- Adding Eqs. B18 and B20 produces 



the result given in Eq. 35 
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FIG. 1. Periodic orbits in the asymmetric infinite square well. The Newtonian orbit N bounces 
back and forth between the infinite walls at x = ±a. The non-Newtonian orbits L and R reflect 
when they reach the discontinuity in the potential at x = and are confined to the left and right 
sides of the well, respectively. 



20 



tg. 400 




FIG. 2. Comparison of approximations with numerically computed eigenvalues. Open squares 



show eigenvalues obtained by numerically solving Eq. 41 Open circles show the perturbation 



theory approximation of Eq. 14 Filled diamonds show the periodic-orbit theory approximation of 



Eq. 37. Parameter values were chosen so that a = 450. 



Necklace 


n L + n R 


r(p) 


xip) 


L 


1 





1 


LRL=(LR)+L 


3 


2 


3 


LRLRL=2(LR)+L 


5 


4 


5 










(i-l)x(LR)+L 


2j-l 


2j-2 


2j-l 



TABLE I. Parameters for periodic orbits with cr(p) = 1 and ni — hr = 1. 
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Necklace 


riL + riR 


r(p) 


x(p) 


R 


1 





2 


RLR = (RL)+R 


3 


2 


4 


RLRLR = 2(RL)+R 


5 


4 


6 










(i-l)x(RL)+R 


2j-l 


2j-2 


2j 



TABLE II. Parameters for periodic orbits with u{p) = 1 and 
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